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Abstract 

Based on the known implicit solution for nonlinear plasma 
waves, an explicit solution was obtained in the form of decomposi- 
tion into harmonics. The solution obtained exhibits a mechanism 
for steepening of nonlinear plasma wave as a result of increasing 
contribution of harmonics and may be used in theoretical studies 
of processes involving nonlinear plasma waves. 
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In the present work an explicit solution for one-dimensional nonlinear 
plasma waves was obtained. Below, prior to the discussion of the solution 
per se, we shall give a formal mathematical description of the method for 
obtaining explicit solutions in case when the implicit solutions are known. 

1) In many nonlinear problems of physics the solutions are sought in 
expressions of the following form 



where f(y) is some nonlinear function. The variable x may be e.g., the coor- 
dinate or the time, while as y the electric field or magnetic field strengths, the 
velocity of the continuum etc. may serve. The expressions of the type (1) are 
usually termed as implicitly given with respect to the quantity y. This equa- 
tion may describe, for instance, some nonlinear oscillation phenomena, the 
results of experimental observations, the distribution of electromagnetic field 
in a cavity or the time dependence of electric field strength in the receiving 
aerial. The implicit solutions of the type (1) are usually highly inconve- 
nient and the knowledge of an explicit solution of the type y — y (x) is much 
preferable. It is frequently required to have an explicit solution in the form of 
decomposition into harmonics or of a spectrum (e.g., as a series or the Fourier 
integral), a series expansion in some set of functions etc. The method that 
is described below will permit the reconstruction from the relation (1) of 
explicit representations for some function g(y) (in particular g — y) in the 
form of series expansion or an integral. 

Let some function q(x) be given. The transform of this function is called 
the quantity (see, e.g., [1]) 



where L x is an operator that usually stands for the integration over the x 
variable, and the function K + (a, x) is termed as the kernel of transformation. 
If Q(a) is known, then one can determine the function q(x) by means of an 
inverse transformation 



x = f(y) 




Q(a) = L x [K + (a,x)q(x)], 



(2) 




(3) 
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In (3) L a is the operator of inverse transformation implying the summation 
or integration over the parameter a, K_(a,x) is the kernel of inverse trans- 
formation. The expressions (2) and (3) describe the Fourier, Hankel, Laplace, 
Mellin transforms etc. Suppose that (1) is given and it is required to find an 
explicit expression for some given function g(y), i.e., to find g = g(x). We 
have from (1): 

L a [K4a,x)R(a)] = L a [K.(a, f)R(a)} (4) 

Here R(a) is not known. The right hand side of the expression (4) is a 
transform of the type (3) for some function of y. We shall require that this 
function be the one of interest to us, g(y). We can then determine the un- 
known function R(a). Thus, 



g(y) = L a [K-(a,f)R(a)], (5) 

R(a)=L f [K + (a,f(y))g(y)]. (6) 

Since the left hand side of (4) is also equal to g(y), then taking into account 
(6) we obtain the required explicit representation for g(y) as a function of x: 

g(x) = L a [K_(a,x)L f [K + (a,f(y))g(y)}. (7) 

So, if we have the implicit expression of the form (1), we can find the explicit 
representation g{x) for some function g(y) by means of the formula (7). 
Although the method under consideration in this or that form is known in 
mathematics, e.g., in the theory of Bessel functions [2], however, it is rarely 
applied in physics. The approach similar to the one discussed above may be 
of help also at calculations of intricate sums over Bessel functions [3]. 

2) Now, we shall determine the explicit solution for the nonlinear one- 
dimensional plasma waves in cold plasma. These waves are described by a 
set of equations consisting of the equation of motion and continuity equation 
for plasma electrons as well as of the Poisson equation: 



— - = -4tt |e| (ra - n e ). (10) 

where t> e and n e are the velocity and density of electrons, n is the equilibrium 
value of density, f - the electric potential that is related with the electric field 
strength by means of the formula E = —dip/dZ, e and m e being the electron 
charge and mass respectively. For steady nonlinear waves propagating with 
the phase velocity v p h, one can obtain from (8)-(10) the following equation 
for the potential [4]: 



^k.^ + 4nn lei 



1 _ V P h 



(2\e\<p/m e ) 1 / 2 \ °' (U) 

where r = uj pe (t — Z/v p h), 0J pe = (Airn e 2 /m e ) l l 2 is the electron plasma 
frequency. Multiplying (11) by dip/dr one can find the first integral in the 
form (dip/dr) 2 / /2 + U(ip) = const. The integration of the latter expression 
gives an implicit solution of the problem [4] (see also [5]): 

r = f(y) = - arccosy - A(l - y 2 ) 1 ' 2 . (12) 

In (12) -1 ^ y = (^ 2 - I)/ A ^ 1, (1 - A) 2 < $ = 2\e\^/{m e v 2 ph ) ^ 
(1 + A) 2 , $ is the dimensionless electric potential , A = v m /v p h ^ 1, where 
v m is the maximum velocity of electrons in the wave. We try to obtain the 
explicit expression for the dimensionless potential <3>(r) as an series expansion 
into harmonics , i.e., in the form of Fourier series (note that the procedure 
to be described below may be pursued also for another quantity, e.g., for 
$ x / 2 ). Since according to (12) $(r) is an even function, we shall seek the 
explicit expression in the form of Fourier series that comprises the cosines 
only. Then, pursuant to the above method, we have from (12), 

oo oo 

a n cos(nr) =^ a n cos(nf) = $. (13) 

n=0 n=0 

Hence, according to the Fourier transformation 
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a n = (2/n) / $cos(n/)d/= (14) 



o 



(2/tt) J (l + Ay) 3 (l-y 2 ) 1/2 cos{n{&rccosy + A(l-y 2 ) 1/2 ]}dy 

-l 

Making the substitution y = cos\l/ we have from (14) 

a n = -(2/tt) y (1 + Acos^) 3 cos[n(^ + Asin^)]rf^. (15) 



o 



Using the well known relations for Bessel functions of integer order (see, e.g., 
[1]) 

u-x) = (-irux), 



TV 

J n (x) = (1/tt) J cos(n^ - xsm^)d^, 
o 

as well as the known recurrent expressions for the Bessel and trigonometrical 
functions, we obtain from (15) 

a n = -(-l) n (A/n 2 )J n (nA). (16) 

Then, according to (13), the required expansion into harmonics for the di- 
mensionless potential will take the form: 

$(r) = -4 jr, (-1)™^M1 cos(nr). (17) 
And for the strength of the electric field we have an explicit expression 



cv x dtp m e u pe v ph d§ m e u pe v ph J n (nA) . 

E(r) = ——- = ; , — — = r-r- 1 — • 2 > —1 sin nr . 

dZ 2 |e| dr e ^ V ; n V ; 

n=l 

(18) 



The nonlinear plasma waves excited by bunches of charged particles were 
formerly studied by simulation of the initial equation (8) - (10) [6]. The plots 
for the potential and strength of electric field in the plasma wave obtained by 

(17) and (18), and those obtained by means of computer simulation of nonlin- 
ear plasma wave excitation [6] practically coincided in case when the ampli- 
tudes of waves were equal . For weak plasma waves, when A = v m /v p h <C 1, 
the series expansions in (17) and (18) are rapidly converging. Hence, in the 
linear case we arrive at the well known result: ip ~ vVax cos r ! E ~ E max sin r. 
The contribution of the harmonics increases with the amplitude and leads to 
characteristic steepening of the wave and eventually to the wavebreaking at 
-4 = 1, E max = m e u pe v ph /\e\. 

Thus, the nonlinear variation of the shape of plasma wave is connected 
with the increase in contribution of harmonics. In the theory of processes 
involving weakly nonlinear plasma waves one can confine to the consideration 
of a few (say, the first two or three) terms in the series expansions (17) and 

(18) . 
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